Multipartite entanglement detection from correlation tensors 
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We introduce a general framework for detecting genuine multipartite entanglement and non full- 
separability in multipartite quantum systems of arbitrary dimensions based on correlation tensors. 
Regarding genuine multipartite entanglement our conditions are comparable to previous approaches 
in the case of qubits while they show particular strength in the relatively unexplored case of higher 
dimensional systems. In the case of non-full separability our conditions prove to be advantageous 
in situations where more than two-body correlations are relevant, where most previous conditions 
turned out to be weak. Moreover, they allow for the detection of fully bound entangled states. 
Finally, we also discuss experimentally-friendly ways of implementing our conditions, which are 
based on directly measurable quantities. 
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I. INTRODUCTION 

In many-body quantum physics entanglement con- 
stitutes a fundamental feature. Complex systems with 
multipartite quantum correlations can be exploited to 
enable numerous tasks in quantum information process- 
ing. The multipartite entanglement in these systems 
enables quantum computation (e. g. multi-party 
cryptography (e. g. @t1|) and the implementation of 
various other quantum algorithms (e. g. (H). Apart 
from these possible applications in modern quantum 
technologies it has become apparent that multipartite 
entanglement also plays a fundamental role in the physics 
of complex systems. While the involvement in quantum 
phase transitions (e. g. [f|) and ionization procedures 
(e. g. Q) seems clear, the recently suggested role in 
biological systems is still subject of debate (e. g. UHllj]). 

Therefore, to decide if a state is entangled or not is 
a fundamental problem in quantum information theory 
[lU, Though a simple mathematical characteriza- 

tion is elusive (the problem has in fact been proved to 
be NP hard [1J]), several works have put up sufficient 
conditions to identify a multipartite state as entangled 
These conditions are not only helpful for en- 
tanglement detection but also they provide more phys- 
ical insight into this phenomenon. Contrary to the bi- 
partite case, there exist different classes of multipartite 
entangled states. Genuine multipartite entanglement is 
of particular interest since it involves entanglement be- 
tween all the subsystems. Recently there has been a lot 
of progress concerning its detection, mostly using linear 
and non-linear entanglement witnesses (l9l - [27l | and Bell- 
like inequalities [28L [23|. However, with a few exceptions 
(see e. g. [13]), the approaches taken in each particular 
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case only allow to detect either entanglement or genuine 
multipartite entanglement. Moreover, most of them are 
limited to qubit systems. 

In this paper we develop a general framework which 
allows using the same piece of information to detect both 
entanglement and genuinely multipartite entanglement 
for multipartite states of arbitrary dimensions. Our main 
tool will be correlation tensors which are built from the 
expectation values of a local operator basis. Our motiva- 
tion stems from different facts. First, it has been shown 
that all information about the entanglement properties 
of a system is encoded in the correlation tensors [30j and 
these mathematical tools have already been proven use- 
ful for the detection of entanglement in the bipartite case 
(3lj . In [l6[ a first step towards the extension of these 
ideas to the detection of multipartite entanglement has 
been taken (see also [32l | for a correlation-tensor approach 
to multipartite entanglement detection). However, here 
we will show that this allows for a much more general 
formalism (in which the criterion of [HI is a particular 
case and that of [32| is strictly weaker), which, further- 
more, enables to identify different classes of multipartite 
entanglement. Our conditions are expressed through sim- 
ple mathematical inequalities. Contrary to entanglement 
witnesses, which are designed for a particular class of 
states, violations of these inequalities signal genuine mul- 
tipartite entanglement or non-full separability for general 
states. Moreover, since the entries of the correlation ten- 
sors are directly related to measurable quantities, we will 
discuss how our approach can be adapted to optimize 
the experimental effort. Last, many conditions for mul- 
tipartite entanglement such as spin squeezing inequali- 
ties [13] , covariance matrice s Hp , entanglement witnesses 
based on structure factors [21] or two-particle Hamilto- 
nians 33] or magnetic susceptibility measurements [33 ] 
rely only on two-body correlations. It has been shown in 
[l8j that this limits their ability to detect entanglement 
as there are important classes of states like graph states 
which have the same two-particle reduced states as sepa- 
rable states. Hence, their entanglement cannot be reveled 
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by just looking at two-point correlations. On the other 
hand, correlation tensors take into account all m-body 
correlations. This suggests (and we will later see) that 
correlation tensors may overcome the limitations of the 
previous criteria. 

II. PRELIMINARIES 

Before we proceed to derive our main results let us 
briefly review the definitions of multipartite entanglement 
and correlation tensors. We consider n-partite quantum 
states p acting on the Hilbert space H = Hi ® • • • ® H n of 
dimension D = d\ ■ ■ ■ d n . If a pure state € H can be 
written as a tensor product of states for every subsystem, 
i. e. 

|*)(*| = |^><V>il®---®IV»»XM (i) 

then the state is said to be fully separable. Consequently, 
fully separable mixed states are convex combinations of 
fully separable pure states. These states contain no en- 
tanglement at all. On the other hand, any n-partite pure 
state that can be written as a tensor product 

\*)(*\ = \* A )(* A \®\* A )(* A \ (2) 

with respect to some bipartition AA (A denoting some 
subset of subsystems and A its complement) is called 
biseparable. These states might contain some entangle- 
ment (as \^a) and/or \^ A ) might not be separable) but 
they are not completely entangled. States that are not 
biseparable with respect to any partition are then said 
to be genuinely multipartite entangled. The generaliza- 
tion to mixed states is straightforward. Any mixed state 
that can be decomposed into a convex sum of biseparable 
pure states is called biseparable. Consequently, any non- 
biseparable mixed state is called genuinely multipartite 
entangled. Due to the fact that the bipartitions might 
differ for every element of the biseparable decomposition 
it is an intricate task to find out whether such a decom- 
position is possible. 

Let {Ap-*}^ 1 denote the generators of SU(dj) and let 

Xq ^ = Idj , which altogether constitute an orthogonal ba- 
sis of the real Hilbert-Schmidt space of Hermitian oper- 
ators acting on Hj (i. e. with inner product (A, B) = 

Tr(AB)). Thus, so is {®J =1 {A^}} for the operators 
acting on H and, hence, p is completely characterized 
by the expectation values (A^ ® • •• ® A^) := % 1 ...i„ 
where ij — 0, 1, . . . ,d? — 1, which gives rise to the so- 
called (multipartite) Bloch representation or (multipar- 
tite) Fano form of density operators |52j. We will de- 
compose the tensor Ti x ...i n into the m-body correlation 
tensors Tj , TjJ k \ etc, which are tensors of order m in- 
dicated by the number of labels in the superscript. All 
the indices not labeled in the superscript are fixed to be 
zero while the other indices take every possible value but 



zero (i. e. the identity is not taken into account). For in- 
stance, the 1-body correlation tensor for particle 1, given 
by = 7iiO - o with ii 7^ 0, completely characterizes 
the reduced state p\ and the 2-body correlation tensor 
for subsystems 1 and 2, T>J 2 ' = Ti^o-.-o ^ 0), 

together with the 1-body correlation tensors of 1 and 2 
characterizes pi2 and so on. For the n-body correlation 
tensor, which we shall also call full correlation tensor, 
we will drop the superscripts to ease the notation, i. e. 
= Ti^-ir, (ij ^ Vj). 
Given two tensors T% x ...i n and Sj 1 .,.j m , their outer prod- 
uct o is the (n + m)th order tensor (To S)ii—i n ji—j m ~ 
Ti 1 ...i n Sj 1 ...j m . If some tensor can be written as the outer 
product of two other tensors, say % 1 ...i n = IZ^^^S^...^, 
we will say that the tensor factorizes in the correspond- 
ing splitting ({1,2,3}, {4, 5, . . . n} in this case). If a ten- 
sor cannot be written as the outer product of any two 
lower order tensors we will say that the tensor does not 
factorize. 

It has been shown in (3lj | that a bipartite pure state 
is separable if and only if (iff) t££ } = T^T^ . Ac- 
cordingly, a multipartite pure state is biseparable with 

respect to the partition AA iff T ( t A) = X {A) T, {A) . Thus, 
we have the following characterization of biseparable pure 
states: 

Fact 1 A pure state is biseparable iff there exists some 
partition of the subsystems AA for which all the m-body 
correlation tensors involving k particles from A and m—k 
from A (k =^ 0,m) factorize into the corresponding k- 
body correlation tensor of the k particles from A and the 
(m — k) -body correlation tensor of the m — k particles from 
A. 

This leads to a simple sufficient condition for genuinely 
multipartite entangled pure states: 

Corollary 1 If some m-body correlation tensor of a pure 
state cannot be factorized into meaningful lower order 
correlation tensors, then the state contains genuine mul- 
tipartite entanglement. 

Analogously, this can be extended to non fully separa- 
ble states (see also [la|): 

Corollary 2 If some m-body correlation tensor of a pure 
state cannot be fully factorized into meaningful 1-body 
correlation tensors, then the state is not fully separable. 

We stress that the factorization must be possible into 
meaningful correlation tensors. This a consequence of the 
fact that the Bloch representation holds for Hermitian 
operators and not only for density operators, which are 
furthermore positive semidefinite. Hence, not all values 
of 7i 1 ...i n give rise to a density matrix, i. e. are mean- 
ingful. To characterize this subset is a quite involved 
problem (see e. g. [35^). However, there exist several con- 
ditions the set of meaningful correlations should fulfill. 
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For instance, it will be useful later on that for 1-body 
correlation tensors it must hold that 



||T«|| < 



'2(di-l) 



(3) 



with equality iff the state is pure and where 1 1 • 1 1 is the 
standard Euclidean norm for vectors. This expresses the 
fact that Trp 2 < 1. This condition holds for the follow- 
ing choice of normalization for the generators of SU(dj): 
Tr(A m A n ) = 2S mn (of course m,n/0 since for the iden- 



tity we have Tr(Ao An 



\q J ') = dj). We will follow this con- 
vention throughout the paper, with which for qubits the 
generators correspond to the standard Pauli matrices. 

Following 0, [3l[ , the main idea behind this paper is 
to express the factorizability of some tensor into lower or- 
der meaningful tensors as an upper bound on some con- 
vex function. Convexity will then imply that this bound 
must hold as well for biseparable (fully separable) mixed 
states and, hence, a violation of this bound will signal the 
presence of genuine multipartite entanglement (non full 
separability) for general quantum states. It seems that 
some tensor norm is the best choice of convex function 
since the norm of meaningful correlation tensors is upper 
bounded as we have just seen for 1-body correlation ten- 
sors. Notice that convexity in this case is guaranteed by 
the triangle inequality. Physical intuition suggests that 
full correlation tensors should be the first ones to check 
and usually we will restrict ourselves to them. 

The rest of the paper is organized as follows. In Sees. 
IIIII and IIVI we provide two different approaches that lead 
to conditions for the identification of genuine multipartite 
entanglement. In Sec.|V]we show that similar techniques 
can be used to obtain conditions for non full separability. 
Section I VII is devoted to some mathematical properties 
of our conditions which are related to their experimental 
implementation. Final conclusions are drawn in Sec. lVIll 



III. GENUINE MULTIPARTITE 
ENTANGLEMENT CONDITIONS BASED ON 
THE STANDARD TENSOR NORM 

The standard tensor norm is defined as the natural gen- 
eralization of the Euclidean vector norm to higher order 
tensors (recall that we will always deal with real tensors) , 
i. e. 



IT,;, 



E 



(4) 



the 1-body correlation tensors. Now, combining this con- 
ditions with Trpfj < 1 will yield an upper bound for the 
2-body correlation tensors (see e. g. [3g|). This procedure 
can be recursively applied to upper bound the standard 
norm of all meaningful rn-body correlation tensors. For 
instance, this gives 



(5) 



with equality iff the state pjk is a maximally entangled 
state [3q . 

To illustrate this, let us start by considering a tripartite 
pure state with subsystems of equal dimension dj = d V j. 
Then, full separability implies 



|T (1) ||||T (2) ||||T (3) | 



2(d- 1) 



3/2 



. (6) 



and biseparability between any two subsystems and the 
other yields 



IT,;, 



li r «,l,, r (*.o,,<^EJM^. (7 ) 



Since the last condition is more restrictive and using con- 
vexity we then have 

Theorem 1 If for an arbitrary (pure or mixed) tripartite 
state it holds that 



IT,, 



> 



8(d- l)(ri 2 - 1) 



(8) 



then the state is genuinely multipartite entangled. 



This seems to be a very good choice for our purposes 
since this norm is multiplicative under outer products, i. 
e. \\T o S\\ = \\T\\\\S\\VT,S. Hence, we will just need 
to upper bound the standard norm of the rn-body corre- 
lation tensors. This turns out to be quite easy. As we 
mentioned above, the condition that Trp 2 < 1 must hold 
translates into an upper bound for the standard norm of tations \\T^ 



This simple mathematical idea is already strong 
enough to detect paradigmatic cases of genuine multipar- 
tite entanglement. If we consider © for three qubits, this 
gives the bound -s/3 ~ 1.73 while the GHZ and W states 
have respectively ||T , j 1 j 2 j 3 || = 2 and IIT^^H ~ 1.92. 
Therefore, their genuine multipartite entanglement is suc- 
cessfully identified, leading to (modest) white noise toler- 
ances of pghz < 0.13 and pw ^ 0.10 fUJ. Furthermore, 
the power of this condition increases with the subsystem 
dimension improving remarkably on [23| . In Figure [1] we 
plot the detection power of Eq. ([8]) for dimensions 4, 5 
and 6. As the other criteria for genuine multipartite en- 
tanglement in high dimensional systems are still based on 
qubit subsystems of the high dimensional Hilbert space it 
is perhaps not surprising that our criteria, exploiting all 
degrees of freedom, quickly outperform them with grow- 
ing dimensionality of the system. 

As discussed above the extension of this condition to 
states with more subsystems or different subsystem di- 
mensions is straightforward. However, for four qubits we 
have for both the GHZ and the Dicke state with two exci- 



= 3, which is precisely the same value 
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Figure 1: (Color online) Here the parameter regions for which the state p — apGHZ(d) + Ppw(d) + 1 ~ id-2 (Xlf=o 

i, i + 1| + |i + 1, i + 1, + 1, i + 1, i\) for (a)d = 4, (b)d = 5 and (c)d — 6 exhibits genuine multipartite entanglement are 
identified. The generalized GHZ and W states for d-dimensional systems are defined as PGHZ(d) = \GH Z(d)){GH Z(d)\ with 
\GHZ(d)) := ^Etol*,i,i> and PwW = |W(d)>(W(d)| with \W(d)) = -^i— £fr o 2 (|i,M + 1> + I*.* + + I* + W))- 

The (red) region labeled II, uses criterion II from Ref. [53] optimized numerically over all local unitary representations of the 
density matrix. The (yellow) region labeled III, uses criterion III from Ref. [23|] optimized numerically over all local unitary 
representations of the density matrix. The numerical optimization was performed using the composite parametrization from 
Ref. [3^ | . The (blue) region labeled C, shows the states detected to be genuinely multipartite entangled using Eq. (JHJ. 



of a tensor product of two maximally entangled bipartite then 
states (||Tj 1 j 2 j 3 j 4 || = v3"\/3). On the analogy of ([8]), one 
might hope to improve for larger d; nevertheless, in the 
next section we will present a different and more powerful 
approach. 



It is worth mentioning that ||Tjj 



[2(d- 



has been shown to be an entanglement monotone [37 
Our results show that a high value of this measure can 
not only imply some entanglement but even genuine mul- 
tipartite entanglement. 



i 

k 

\\A\\ k =Y J ^ 



(9) 



IV. GENUINE MULTIPARTITE 
ENTANGLEMENT CONDITIONS BASED ON 
NORMS OF MATRICIZATIONS OF TENSORS 



As an alternative to the previous section one can seek 
for other norms. Unfortunately, to our knowledge, the 
standard norm is the only norm which is multiplicative 
under outer products, a property which is very convenient 
for the mathematical simplicity of our derivations. Nev- 
ertheless, it turns out that considering matricizations of 
tensors (i. e. particular rearrangements of the tensor val- 
ues to form a matrix) [39| and the usage of matrix norms 
on these matricizations will lead to interesting and more 
powerful results. The matrix norms we will be dealing 
with are the Frobenius or Hilbert-Schmidt norm (which 
is the standard tensor norm on a matrix), the trace norm 
and the Ky Fan k norms [40|. That is, let A G M mx ", 



where {(Ji\ (i — 1, . . . , min(m, n)) denote the singular 
values of the matrix, which are arranged, as usual, in 
non-increasing order. Notice that the last Ky Fan norm 
is the trace norm, i. e. || • || m in(ro,n) = II ' lltr- 



We will define matricizations in the following way: 
(non) underlined indices are joined together in lexico- 
graphical order to give rise to the row (column) indices. 
Let A be the subset of underlined indices, then we will call 
that a A, A matricization. For example, let ij = 1, . . . , rij 



■5 



Vj, then 



/ Tnik 

T\12k 



T\21k 
T\22k 



21n 3 k 



Tn 



(10) 



kJ 



where T- 



xyzk 



yzni) is a row vector (fc 



1,2,..., 71,4), is a 13,24 matricization. In Dirac notation 
we would have 



|«1»3)(*2«4| 



(11) 



This way of matricizing is a generalization of the concept 
of matrix unfolding or mode-n matricization which is of- 
ten used in multilinear algebra 391], which corresponds 
to matricizations of one index giving rise to the row col- 
umn and the rest to the column vectors. The matri- 
cizations we have defined are quite convenient for the 
problem at hand because they have a well defined struc- 
ture under outer product of tensors. Of course, if all the 
indices of a tensor are joined together the tensor is vec- 
torized, T^i^ = uec(T), while it is straightforward to 
check that T^W^ = (T ili2 ) (g) (Wj 3 i 4 ). Concatenating 
these rules a matrix form for any matricization of more 
involved outer products of tensors can be readily found. 
For instance, 

= (Ti 2 ,ni 3 ) ® (^ 4 ) T ® vec{S l5l6 ) ® (W lTlg f . (12) 

This is the kind of structure we need because the norms 
we are going to use are either multiplicative ( 1 1 • 1 1 , 1 1 • \ \t r 
and I Mil) or submultiplicative (| | • | |fe) under tensor prod- 
ucts [54]]. So, analogously to the previous section, we just 
need to upper bound these quantities for meaningful cor- 
relation tensors to obtain conditions for genuine multi- 
partite entanglement. For the sake of simplicity we will 
consider multiqubit systems in the following subsections. 



A. Three qubits 

According to the above discussion the only thing left to 
be able to derive genuine multipartite entanglement con- 
ditions is to obtain upper bounds for the matrix norms of 
the correlation tensors similarly as we did with the stan- 
dard norm in the previous section. We have the following 



Lemma 1 The 2-body correlation tensor of two qubits 
satisfies 



\T^\\k<k Vfc, 



(13) 



with equality iff the two qubits are in a maximally entan- 
gled state. 

Proof. We will use the local unitary invariance of the 
norms of any matricization of the correlation tensors (see 
Sec.rVTlbelow). Notice then that, due to 50(3) ~ 517(2), 
the 2-body correlation tensor (i. e. correlation matrix) can 
be brought into diagonal form by choosing properly lo- 
cal unitaries in the two subsystems (see e. g. [4l|). Since 
the entries of this matrix are expectation values of ob- 
servables with eigenvalues 1 or — 1 we have that \Tu\ < 1 
Vz. It can be readily checked that this bound is attained 
by the maximally entangled state (and only by the maxi- 
mally entangled state because this value of the trace norm 
of the correlation matrix implies the maximal possible 
amount of entanglement, as measured, for instance, by 
the concurrence [42[ ). □ 

Lemma 2 If a pure 3-qubit state is biseparable, then it 
holds that: 



(i) If the state is fully separable 



IT, 



tj'iih, 



< 1 Vfc. 



(14) 



(ii) If the state contains no entanglement across j\lm 



I^VfiltmlU < ^3 Vfc. (15) 



(iii) If the state contains some entanglement across j\lm 



m 



ijii i m 1 1 k 



, < fc Vfc. 



(16) 



Proof. We will use repeatedly the upper bounds (JSJ) 
and (HI 



(i) 



m 



\ k = llT^Tf T^IU = ||(T«) • (T« ® T<"0) r |U 
= ||T w ||||T (i) (g)T (m) || 

= ||T^)||||T«||||T( m )|| = 1. (17) 



(ii) 



= II^^IU = ll(r (i) ) • vec(T^f\\k 
= ||T^||||T^ m >|| < VS. (18) 



(iii) 



\ T i A iiim\\k = \\T^T^\\ k = H(T^)) ® (T^fWk 

< ||T°'' IU||T (m) || < fc. (19) 



6 



□ Lemma 3 The 3-body correlation tensor of three qubits 



From Lemma [5] we read 3 sufficient conditions for gen- 
uine multipartite entanglement, namely that the norm of 
any of the three possible matricizations of the full cor- 
relation tensor is greater than -s/3, 2 and 3 for || • ||i, 
|| • ||2 and || • || tr respectively. To illustrate the power 
of these conditions consider that the singular values of 
these matricizations are {1.414,1.414,0} for the GHZ 
state and {1.374,0.943,0.943} for the W state. Hence, 
the last two conditions can detect genuine multipartite 
entanglement. Notice that these states are symmetric, so 
all matricizations of the full correlation tensor are equal; 
however, for general states our ability to detect a state 
as genuinely multipartite entangled might depend on the 
choice of matricization. To avoid this and to obtain a 
stronger condition which takes into account a combina- 
tion of the bounds of Lemma [2] rather than just picking 
one of them we introduce the average matricization norm 



l|M(T ili2i3 )|| - 
leads to 



, ||)/3, which 



Theorem 2 If for a 3-qubit state it holds that 
\\M(T ili2i3 )\\k>^P^, 



(20) 



then the state contains genuine multipartite entangle- 
ment. 

Proof. Simply use that any biseparable state can be writ- 
ten as p bs = Y.k PkPu ® Pt + 1kPi 3 ®p\ + r k p\ 3 ® p\ and 
combine properly the bounds ([TBI and (1TB)) . □ 
Thus, Theorem [5] allows to detect genuine multipar- 
tite entanglement in mixtures of the GHZ and W states 
with white noise for noise levels of pghz ^ 0.324 and 
Pw 0.209. Notice that it is known that these states are 
genuinely multipartite entangled iff pghz 0.571 and 
p w < 0.521 [23. 



B. Four qubits 

Now, similarly to previous sections, we need upper 
bounds to the norms of the matricizations of the 3-body 
correlation tensor. However, it is not clear which states 
should attain the maximum values of these norms in op- 
position to the 2-body case, where the maximally entan- 
gled state, as intuition would suggest, does the job. More- 
over, numerics indicate that max^ ||TLiii m ||tr — 3.272 
for a state which, although close to the W state, has 
no simple mathematical structure. This indicates that 
devising a systematic procedure to find the maximum 
value of these norms similarly as we did with the stan- 
dard norm in Sec. Mil might be very hard. Nevertheless, 
it turns out that we can use this procedure to obtain rea- 
sonable estimates by using the equivalence of the norms: 

ll-IU< v^lHI. 



satisfies || || fc <2y/kVk. 



Proof. The result follows from \\T t 



ijkim I 



this we proceed as in Sec. [TTTJ The fact that Ttp% m = 1 
translates to 



< 2. To see 

jlm 



J2 l|r (s) || 2 + ^||r^|| 2 -H|r w < m) || 2 = 7. (21) 



s—j,l,m 



s<q 



The minimum possible values of the norms of the lowest 
order correlation tensors are | \ \ — and | |T( s ' g ) 1 1 = 1 
Vs, q (since the 1-qubit reduced density matrices can be 
maximally mixed but the highest mixing allowed by the 
2-qubit reduced density matrices is them being equal to 
the identity in a two-dimensional subspace). Therefore, 
||TW< m )|| < 2, which is attained by the GHZ state. 
Finally, notice that the Hilbcrt-Schmidt norm of any 
matricization of a tensor equals its standard norm as 
tensor. □ 
Notice that Lemma [3] provides accurate estimates as 
the trace norm bound 2-\/3 ~ 3.464 is quite close to 
the numerical maximum given above while the Ky Fan 
2 norm bound is actually sharp since it is attained by the 
GHZ state. 

Now, we can proceed as in Lemma [5] to upper bound 
\\Ti j i l i m i 3 \\k- Then, for pure biseparable states, one ob- 
tains the bounds 2\fk (for k < 3 and 2\/3 otherwise) 
if Ti. 



L ijilt m is — ± i j ± iii m i s UL — ijhi m i s \ 1 ' 

e. the state is biseparable in one subsystem versus the 
other three and we consider the two possibilities that the 
two indices giving rise to the row of the matricization ei- 
ther belong to unentangled or entangled particles), 3 if 

Ti j i l i m i a = T^' k T im l s and k if T ijilirl 



or Ts„ 



rn(j,l,m)rp(s) 

ijhi„ J 



_ rp(j,m)rp(l,s) 

ijim iii s ' 

One could also consider 1 vs. 3 matricizations of the full 
correlation tensor; however, one finds that | |2j.,-jj m j s | < 



y/3k for T ijilimis = T^T^f> (i. e. any 2 vs. 2 bisepa- 
rable state), which turns out to be a weak condition and, 
then, it is better not to take these matricizations into 
account. Combining all the above bounds, defining the 
2 vs. 2 average matricization norm ||M22(T , j li2i3 j 4 ) M — 



\Ti, 



+ \\Ti,i 



+ I |21 lis , i3 i 4 1|)/3 and proceeding 



as in Theorem |2| we have 

Theorem 3 // for a J^-qubit state one of the following 
inequalities holds 



\M 22 {T ni 



k > 



( 2y/k 1 < k < 3 

\ l + 2fc/3 4 < k < 9 



, (22) 



then the state contains genuine multipartite entangle- 
ment. 

With this, genuine multipartite entanglement is de- 
tected in the GHZ state with a white noise tolerance of 
Pghz ^ 0.307 and for the Dicke states of 1 and 2 excita- 
tions we have respectively pr> x < 0.018 and pd 2 0.328. 
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From (271 ] we know that there is genuine multipartite en- 
tanglement iff pghz i$ 0.533 and if po 2 < 0.539. 

These examples indicate that the matricization ap- 
proach is more powerful than that of the standard norm 
and that it can detect different classes of entangled states. 
Interestingly, the detection capability of Theorems [5] and 
[3] is already comparable to [23J, [26[ for qubits as shown in 
Figure [2j 
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Figure 2: (Color online) Here the parameter regions for 
which the state p = ap GHZ{2) + f3p D ( 2 ) 



hibits genuine multipartite entanglement are identified. The 
GHZ and two-excitation Dicke state for four qubit sys- 
tems are defined as p G HZ(2) = \GHZ{2)){GHZ{2)\ with 
\GHZ(2)) := ^ELoK-M) and p D{2) = \D 4 2 }{D 4 2 \ with 
\D 4 ) = ^=(j0011) + |0101) + |1001) + |1010) + |1100) + j0110}). 
The (red) region labeled II, uses criterion II from Ref. [xj op- 
timized numerically over all local unitary representations of 
the density matrix. The (yellow) region labeled I 2 , uses cri- 
terion I 2 from Ref. [26l ] optimized numerically over all local 
unitary representations of the density matrix. The numerical 
optimization was performed using the composite paramctriza- 
tion from Ref. [3^|. The (blue) region labeled C, shows the 
region detected to be genuinely multipartite entangled using 
theorem [3] 



Using the matricization approach we have thus con- 
structed versatile criteria, detecting genuine multipar- 
tite entanglement in a broad variety of cases. All fa- 
mous examples of four qubit multipartite entangled states 
are detected (GHZ-, W-, Dicke- and Singlet-state) using 
the same criterion, without any optimization involved as 
the norms of any matricization of a correlation tensor 
is invariant under local unitary transformations on the 
density matrix (see Sec. [VTJ). Although for some specific 
states optimizing over all possible witnesses can yield a 
higher noise resistance in some cases, a comparable result 
is achieved in a computationally far more efficient way. 
As shown in Figure [3J there even exist states that were 



not detected to be genuinely multipartite entangled with 
any of the optimized criteria so far. 



V. DETECTION OF NON-FULLY SEPARABLE 
STATES 

As mentioned above, we can also use correlation tensors 
to discriminate states containing some form of entangle- 
ment and fully separable states. This has already been 
carried out in [lfij], where the authors show that that for 
fully separable states an upper bound on the trace norm 
of the full correlation tensor must hold for any matriciza- 
tion of the form one particle versus the rest (i. e. matrix 
unfoldings). However, remarkably, our picture allows not 
only for a very simple proof of this fact, but, also, for a 
significantly stronger result since we can show that such 
a bound must hold for any possible matricization of the 
correlation tensor. 

Theorem 4 For any fully separable state any matriciza- 
tion of the full correlation tensor must fulfill 



-il---ikik-t-l---in \\tr 




i=i 



i. e. k = 1, . . . , n — 1 and all possible permutations of the 
particles are taken into account. 

Proof. Since the correlation tensor of a fully separable 
pure state must fully factorize into the 1-body correlation 
tensors it is straightforward to note that 



IT, 



il---ikik + l---in \\tr 



, T (1) _ _ T (fe) T (fc+l) _ _ . r (")| 
I j± ik_ ik+l in ' 



= ||(T (1) ® • • • ® T«) • (T( fe+1 > ® • ® T(")) T || tr 

= ||T (1) ® ••• ® T (fc) || \\T^ k+1) ® ••• ® T (n) || 

= Y[\\T^\\ (24) 

3 

Using Eq. the proof is finished. □ 
Notice that similar bounds hold as well for lower order 
correlation tensors. Of course, one could also consider 
other matrix norms but the trace norm yields the most 
powerful condition as they all lead to the same bound. In 
particular, if one considers the standard norm one obtains 
the criterion of [32| . which is then proved to be strictly 
weaker than our Theorem [3] as 1 1 • 1 1 t r > \ \ ■ 1 1 . 

Theorem 2] is then clearly stronger than the criterion 
of [iH | since it contains matrix unfoldings as a particular 
case while the other matricizations can further restrict 
the set of fully separable states leading to a substan- 
tially more powerful detection of entangled states. For in- 
stance, consider the Dicke states of four qubits for which 
the matrix unfolding (i. e. the criterion of |16j ) detects en- 
tanglement in some form up to white noise levels of 0.698 
(1 excitation) and 0.807 (2 excitations), while the 12,34 
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matricization rises these levels to 0.732 and 0.842 respec- 
tively. As we will see in more detail with the examples 
below this seems to be a general feature. 

A very interesting conclusion of our study of multi- 
partite entanglement detection with correlation tensors 
is the fact that the very same piece of information can 
be used to decide both non full separability and gen- 
uine multipartite entanglement. This gives Theorem H] 
an advantage over other conditions for the detection of 
some form of entanglement in multipartite states, since 
used together with Theorems [5] and [3] one can further- 
more discriminate when genuine multipartite entangle- 
ment is present. Moreover, as discussed in the introduc- 
tion, the correlation tensors take into account all m-body 
correlations while many important conditions rely only 
on two-point correlations. This constrains the power of 
these conditions since, for instance, they cannot detect 
the important class of graph states as it has been shown 
in 18]. In particular, this means that, contrary to our 
case, they cannot detect the GHZ state, one of the most 
paradigmatic multipartite entangled states. As we will 
see with some examples this limitation is extendable to 
other classes in which interactions of more than two par- 
ticles are somehow relevant (e. g. graph states cannot be 
nondegenerate ground states of Hamiltonians containing 
at most two body interactions |43|). On the contrary, 
Theorem Q] turns out to be quite efficient in these cases. 

Let us start by considering the 4-qubit 3-body inter- 
action Hamiltonian with transversal magnetic field of 
strength h 



H 1 = j2(-^- 1) ^ ) ^ +1) + ^ ) i 



(25) 



where periodic boundary conditions arc assumed and the 
superscript indicates on which qubit the operation is act- 
ing. For brevity, the tensor products and the identity 
operation are omitted. In Figure [3] we plot the detec- 
tion efficiency of correlation tensors and the optimal spin 
squeezing inequalities (OSSI) [lj] for the thermal states 
of Hi [53, i. e. 



p{kT 7 h) 



exp(-Hi/kT) 
Tr(exp(-if 1 /feT))' 



(26) 



The superiority of our condition is not so surprising as 
it is known that the (nondegenerate) ground state of Hi 
when h — is the cluster state, a graph state. Moreover, 
one can consider a slight variation of H±, namely 



for which the thermal ground state when h = is not 
a graph state (5(|. Nevertheless, correlation tensors re- 
main more powerful than the OSSI as shown in Figure 
IH and, furthermore, a region of genuine multipartite en- 
tanglement can be identified with them. Let us mention 




Figure 3: (Color online) Maximal temperature for which 
entanglement is detected in the thermal states of H\ by: The- 
orem 2J i. e. any matricization of the full correlation tensor, 
(x marks), the criterion of [TjJ, i. e. matrix unfoldings of the 
full correlation tensor (pluses) and OSSI (circles) . 




Figure 4: (Color online) The same as in Figure[3]for H2- The 
area in the down-left corner represents genuine multipartite 
entanglement as detected by Theorem 



that these Hamiltonians are not artificial, they can arise 
as effective interactions in the context of optical lattices 
in a triangular configuration (iij and their statistical me- 
chanical properties are a subject of current research (45} . 

In the same vein we have considered the class of locally 
maximally entanglable (LME) states (4(|, which gener- 
alize graph states. These states can be generated with 
nonlocal gates acting on a product state, these operations 
being generalized phase gates with m-body interactions. 
We have checked Theorem 0] in randomly generated 4- 
qubit LME states and they were always found to be en- 
tangled (57(. Moreover, in most cases they were found to 
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be genuinely multipartite entangled by Theorem |3l On 
the contrary, the OSSI always failed for these states [58[ . 

Last, we have checked that Theorem @] is able to de- 
tect entangled states with a positive partial transposi- 
tion (PPT). This was already known in the bipartite case 
[3ll l36j . However, one may wonder if this is still pos- 
sible in the multipartite case for states which are PPT 
with respect to every possible bipartition. This question 
is answered in the affirmative. We have considered the 
3-qubit Hyllus state [47j which depends on one free pa- 
rameter. For all values of this parameter, despite not 
fully separable, this state is not only PPT with respect 
to every possible bipartition but separable with respect 
to every possible bipartition. It can be readily checked 
that Theorem 2] detects this state as entangled for many 
values of the free parameter. This shows the strength of 
our condition since the entanglement in this state is very 
weak: it is fully bound entangled, i. e. the state cannot 
be distilled even if different parties act together. 



VI. MATHEMATICAL PROPERTIES AND 
EXPERIMENTAL IMPLEMENTATION 

As we have shown the correlation tensor provides a 
powerful tool, enabling the detection of a wide range of 
multipartite entangled states. It is worth pointing out 
that 

Fact 2 All the quantities needed to apply the conditions 
presented here ( Theorems {J\^ are invariant under local 
unitary (LU) transformations on the density matrix. 

More precisely, all the norms considered here (|| • ||, 
|| • ||fe and || • \ \t r ) of any matricization of any correlation 
tensor are invariant under these transformations. This is 
because of the well-known fact that LU operations acting 
on the density matrix correspond to rotations in corre- 
lation space (see e. g. [HI HH, EH). Hence, the trans- 
formation p — > U\ ® • • • ® U n pU\ <g> • • • ® U X amounts 



to Ti 



'l-'fu-l"' 1 ' 



-> (Oi 



•O k )T, 



'r-'fe't+i'" 1 ' 



,(o, 



fc+1 



• • • §5 O n ) for any matricization of a correlation tensor T 
for some real orthogonal matrices {Oi}. The result then 
follows because the aforementioned norms are all unitar- 
ily invariant (40j . 

This is a convenient property which is not shared by 
other conditions for genuine multipartite entanglement 
[22I l23l | . Although one can nevertheless take the effort of 
optimizing over LUs for the application of these criteria 
in practice, LU invariance is a very satisfactory property 
from the theoretical point of view since this is a funda- 
mental property of entanglement. Moreover, this could 
lead to the use of these norms of correlation tensors not 
only as qualitative indicators of entanglement but also as 
quantitative tools (this is already the case for the trace 
norm in the bipartite case 42| and for the standard norm 
in the multipartite case for non-full separability (36l.l37l|). 
Furthermore, LU invariance is also quite convenient from 



the experimental point of view since this implies that 
the local measurement settings of each party need not be 
aligned with the others (59|. 

Regarding experimental implementation, it is also 
worth discussing the number of measurements required to 
use our criteria. Although knowledge of the full correla- 
tion tensor requires less parameters than that of the den- 
sity matrix, the statistical data that needs to be collected 
to reconstruct it allows, however, to do state tomography. 
Nevertheless, partial knowledge of a correlation tensor al- 
lows as well to implement our conditions as the norm of 
a matrix can be estimated from below by knowing just 
some of its entries. This is clear for the standard norm 
from Eq. Using again the equivalence of the norms, 
e- g- || • \ \tr > II ' Hi one can lower bound the other norms. 
In addition to this, if particular entries of a matrix are 
known one can directly obtain lower bounds for the Ky 
Fan norms. For instance, we have that 



U\\tr>J2 



(28) 



and, also that the Ky Fan norm of any principal subma- 
trix is a lower bound for the Ky Fan norm of the full 
matrix. Furthermore, the minimum value of the trace 
norm of a matrix subject to the knowledge of some of 
its entries can be efficiently c omp uted using convex opti- 
mization techniques (see e. g. 49]). Hence, if the numbers 
given by these lower bounds are already large enough to 
violate the inequalities given by Theorems HHH one can 
then conclude with certainty the presence of genuine mul- 
tipartite entanglement or non full separability with con- 
siderably fewer measurements. 

Genuine multipartite entanglement is usually ad- 
dressed in terms of entanglement witnesses. They are 
locally measurable and the number of required mea- 
surements can scale very favorably with the system size 
(i.e. polynomially in Refs. [23l . l2q | and even linearly in 
Ref. |25|). However, this limited number of measurements 
of course also severely limits the number of states that 
are detected by such criteria (e.g. GHZ states using the 
criteria from Refs. [23|, |26(). Nevertheless, if one has the- 
oretical expectations of what the state should look like, 
one can then use a suitable criterion which should be able 
to detect these states. Our presented framework can al- 
low as well for a versatile detection of any state that is 
detected by our criteria using only a very limited number 
of measurements as in many cases a considerable number 
of expectation values (i. e. elements of the correlation 
tensors) are zero. 

• First calculate all correlation tensor elements of the 
theoretically expected state in an experiment 

• Second, only perform the measurements corre- 
sponding to exactly these elements 

• Then lower bound the norm of the correlation ten- 
sor using only this limited amount of elements 
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If this lower bound exceeds the threshold of any of our 
inequalities it is certain that it contains entanglement. 
This implies that in cases where entanglement witnesses 
are applicable (i.e. some prior expectations of the state), 
we can apply our criteria in a just as experimentally fea- 
sible way. E.g . for three qubit GHZ states the criterion 
from Ref. 23] requires seven local measurement settings 
which is exactly the number of correlation tensor elements 
that have to be ascertained. 



VII. CONCLUSIONS 

We have provided a general framework to detect dif- 
ferent classes of multipartite entanglement in systems of 
arbitrary dimension using as a main tool correlation ten- 
sors. In particular, considering several norms on these 
objects, we have shown that different upper bounds can 
be established such that violations signal the presence of 
either non-full separability or genuine multipartite entan- 
glement. Regarding genuine multipartite entanglement, 
we have explicitly worked out the case of tripartite qudit 
(Theorems 1 and 2) and four partite qubit (Theorem 3) 
systems. The approach, however, can be generalized to 
an arbitrary number of subsystems and dimensions in a 
tedious but systematical way. It would be interesting to 
study in the future if this procedure can be rendered more 
straightforward. This is the case for our sufficient con- 
dition for non-full separability (Theorem 4), which has 
a very simple proof for arbitrary number of subsystems 
and dimensions. 

We have demonstrated as well with exemplary cases 
that our approach can improve and complement previous 
comparable criteria in both the non-full separability and 
genuine multipartite entanglement scenarios. Further- 
more, the entries of the correlation tensors are directly 



related to measurable quantities and we have discussed 
how to estimate the relevant norms with fewer measure- 
ments to ease the experimental implementation. Last, 
our norms are all LU invariant, which besides some imple- 
mentation advantages, is a satisfactory property from the 
theoretical point of view. This suggests that they might 
be connected to the quantification of entanglement with 
entanglement measures, thus not only providing quali- 
tative information. This would be particularly interest- 
ing in the genuine multipartite entanglement case, where 
the first quantification steps have been taken in [271 l50j . 
We leave for future research the question of whether our 
norms could provide a (rough) easily computable quan- 
tification of genuine multipartite entanglement. 

Another interesting point for future research is the con- 
nection between the norms of the correlation tensor and 
non-locality. It has been shown in that the set of 
two-qubit states violating the CHSH inequality is char- 
acterized in terms of the singular values of the correla- 
tion matrix, i. e. a state violates CHSH iff the sum of 
the squares of the two largest singular values of T^ 1 ' 2 ) is 
greater than 1. It would be interesting to study if this 
connection can be extended to the multipartite regime 
using the norms and correlation tensors we have intro- 
duced here and to establish analogous characterizations 
for other Bell inequalities. 
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